Contour Dynamics Simulations with a Parallel Hierarchical-Element Method

Many two-dimensional incompressible inviscid vortex flows can be simulated with high efficiency by means of the contour dynamics method. Several applications require the use of a hierarchical-element method (HEM), a modified version of the classical contour dynamics scheme by applying a fast multipole method, in order to accelerate the computations substantially. In this paper it is shown that the acceleration of contour dynamics simulations by means of the HEM can be increased further by parallelising the HEM algorithm. Speed-up, load balance and scalability are parallel performance features which are studied for several test examples. Furthermore, typical simulations are shown, including an application of vortex dynamics near the pole of a rotating sphere. The HEM has been parallelised using OpenMP and tested with up to 16 processors on an Origin 3800 cc-numa computer.